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Abstract 

We deal here with the apphcation of the Nuclear Born Oppenheimer 
(NBO) method to the description of nuclear rotations. As an edifying il- 
lustration, we apply the NBO formalism to study the rotational motion of 
nuclei which are axially-symmetric and even, but whose shells are not closed. 
We focus, in particular, on the derivation of expressions for the rotational 
energy and for the moment of inertia. Additionally, we examine the con- 
nection between the NBO method and the self-consistent cranking (SCC) 
model. Finally, we compare the moment of inertia generated by the NBO 
method with the Thouless-Valantin formula and hence establish a connection 
between the NBO method and the large body of experimental data. 

PACS numbers: 21.60.-n, 21.60.Ev, 21.10.Re 
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1 Introduction 



Since nuclear and molecular rotation-vibration spectra present many striking analo- 
gies, and since the Born-Oppenheimer (BO) approximation[l] of molecular physics 
was shown to be very accurate^ in describing molecular rotations and vibrations[2], it 
will be interesting to explore the possibility of using the BO approximation to describe 
nuclear collective rotations. 

Exploiting the analogy between nuclear and molecular dynamics, Villars introduced 
a microscopic method[3], [4] to describe nuclear collective motion. This method, to 
be called the Nuclear Born-Oppenheimer (NBO) method, was developed along the 
lines of the molecular BO approximation by constructing a factorable trial function 
modeled after the BO ansatz. 

Using an analytically solvable model[5], we have shown that the NBO method is 
very accurate for adiabatic collective motion[6]. Since the NBO method is a quan- 
tum mechanical prescription, we have shown that it offers a suitable framework for 
describing the zero-point fluctuations[7]; we have also shown that the method offers 
an accurate description of small-amplitude collective oscillations[8] and that it yields 
the random phase approximation (RPA) equations[9]. Additionally, we have applied 
the NBO method to study nuclear collective motion[10] and examined its connection 
with the collective model of Bohr[ll]. 

So, having applied the NBO method to the study of small amplitude motion, we 
have yet to apply it to nuclear collective rotations. In this work we want to achieve just 
that aim; namely, we want to apply the general BO formalism outlined in Ref.[10] to 
the description of nuclear rotational states. As an illustration, we will apply the NBO 
formalism to study the rotations of nuclei that are axially-symmetric and even, but 
with non-closed shells. We will focus, in particular, on the derivation of expressions 
for the energy and for the moment of inertia. Additionally, we shall examine the 
connection of the NBO method with the successful self-consistent cranking (SCC) 
model. 

In Sec. 2, we present a brief outline of the NBO formalism and how it applies 
to nuclear nuclear collective motion. We then devote Sec. 3 to the application 
of the NBO method to study the rotational states of axially symmetric nuclei; in 

^Using an elementary solvable model, Moshinksy and Kittel[2] have shown that the BO 

approximation is very accurate for both the molecular energy and wave function: ~ 

and |{V'Bo|'^^exact)P = 1 — ilfX^i where x is equal to the ratio of electronic to nuclear masses 

(i.e., X = rWe/M ~ 10*). 
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particular, we will derive an expression for the rotational energy. In Sec. 4, we 
present a discussion on the connection of the NBO method to the self-consistent 
cranking (SCC) model. 

2 Synopsis of the Application of the NBO Method to 
Nuclear Collective Motion 

To describe nuclear collective motion within the framework of the NBO method, we 
need to introduce a tensor operator Qai3'< that is, to be able to describe collective 

rotations and vibrations of nuclei, we need to introduce a set of operators that are the 
elements of a symmetric cartesian tensor operator Qap- In the rest of this work, we 
shall use Greek subscripts to refer to a space-fixed frame of reference; Latin subscripts 
will be used later to refer to a body-fixed frame. The operator Q^/s is assumed to 
depend on the various nucleonic variables - positions xi, momenta, pi, and spins, Sj. 
In addition, we assume that Q^fS are one-body operators, symmetric, even under time 
reversal, have a continuous eigenvalue spectrum (qajs), and commute with any other 
component Qjs of Q { i.e., [Qa/B, Qjs] = 0)- Let Ka/s be the canonical conjugate 
of Qal3- 

[iKap, Q-ys] = 2 {^a-ySps + SasS/S-y) ■ (2.1) 

The NBO method consists of the following two essential steps[3], [4], [10]: 

• First, we need to construct a suitable representation for the nucleus' Hamilto- 
nian H by decomposing it into a series 

H = Ho + ^Hl apKaP + ^ ^2 aP^sKapK^S + " " " , (2-2) 

a/3 aj3 7<5 

where all coefficient operators Hq, Hia/3, H2ai3,jS, ■ ■ ■ commute with Qa/3- 

• Second, we make use a factorable trial function^ 

= / n dQap{^\S{Qap - Qa(3)Mq))g{q) , (2.3) 

0!</9 

where {x\^{q)) is the intrinsic wave function, and g{q) is the collective ampli- 
tude. 

^We will use x to abbreviate for the set of nucleonic variables - position Xi, momentum pi, 
and spin, Sj. 
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After constructing the Hamiltonian and the wave function, we can calculate the 
mean energy by a simple application of H to IV'): 



\H\ 



a<(i al3 



1 



+ 9 S!S!^2a/3,76|^)W^65'(«) + •••}, 



(2.4) 



aj3 -yS 



where the ka/3 are operators that act on g{q); they obey commutation relations with 
the Qap isomorphic with (2.1) 



[ikaf3, q-ys] = 2 i^aySps + SasSfiy 

The few lowest expressions of Hk are given by 

Ho = H - J2 iQap] Ga/J + ^ EE [[^' 

a/3 a/3 jS 



(2.5) 



a/S 



Hi a/3 = 



H2a/3 jS- 



(2.6) 



7<5 



H, iQa/3] , iQjs] + 



where Gq/? is a one particle operator that acts on |0); it is defined by the action of 
ka/3 on the parameter q in 10) 



1 d 

kap{x\^{q)) = ;rT(l + (5a/3)^ (^l^(^)) = $((?)). 

2l aqa/3 



(2.7) 



We should note that the mean energy expression (2.4) was derived within a space- 
fixed or lab frame. However, in the description of permanently deformed (non spher- 
ical) nuclei, it is more convenient to employ a body-fixed frame of reference. Here, 
we take the axes of the body-fixed frame along the three principal axes of qap which 
are defined by the unit vectors Ba (a = 1,2,3), and specify their orientation with re- 
spect to the space-fixed frame by three Euler angles[23] 9s (i.e., 9,(f,il)): = ea(9s) 
with ia ■ et, = 6ab and ia x eb = Sahc&c where £abc is the antisymmetric tensor 
(^123 = 1 = —^213 etc.). The collective degrees of freedom can be separated into 
rotational and vibrational terms by transforming Qa/3 to the body-fixed frame; that 
is, by means of the principal axes transformation of the tensor operator Qa/3'- 



qa/3 = Y^eaa{0)ei3a{9)qa, 
a=l 



(2. 
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where Caai^ ■ Ca). the a component of the unit vector ea, depends on the three 
Euler angles Bg- In the transformation to the body-fixed frame, we have essentially 
replaced the six collective coordinates g^/? by the three qa% and the three Euler 
angles. The matrices eo.fi obey the orthogonality relations: Yl,a^o-a^fia = and 

We can now introduce rotation operators about the body-fixed axes. The 
Euler angles specifying the orientation of the intrinsic frame need to be viewed as 
dynamical variables; for instance, the unit vector Ca satisfy the commutation rules of 
a vector operator 



C 



ab] ) 



(2.9) 



We can easily verify from (2.9) that these operators obey the commutation relations 



(2.10) 



which differ in sign from the commutation rules of ordinary angular momentum[12] 
because they refer to the moving axes and hence do not have the same commutation 
properties as angular-momentum components along space fixed axes. For instance. 



we have 



—iCs- The space fixed components Cap of £. can be obtained 

by rotation: Ca/s = J2ab 

In conjunction with the replacement of qap by the variables and Og, we seek 
an expression for the operator kap in terms of the jO-ap and a set of three operators 
Pa conjugate to Qa, with \pa, Qb] = iSab- We can verify[10] that ka/s transforms like 
an operator that acts on 6s and qa- 



k 



1 y-v e, 

^ ab 



Qa Qb a 



eaaGpb 

ab 



- (i + dap) 



Qa - Qb 



+ SabPa 



(2.11) 



Additionally, we can ascertain that k^p is Hermitian with regard to the volume element 
Y[a<p'^QaP: which can be shown to transform like: 

3 

j dr = j Yl dQap = / n d'laiqi - q2){q2 - qz){q2, - qi)dO,, (2.12) 



a<p 



where dO, is the usual angular element[23] dil = sin 9 dOdipdtp. 

Using the relations (2.11) and (2.12), we can now express (2.4) and (2.6) in 
the body-fixed frame. For this, note first that under the transformation (2.8) from 
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the Lab frame to the body-fixed system, the quantities g{qa/3), {xiaSia\^{Qa/3)) , and 
Ila</3 ^{Qaf3-Qai3) become f{q \(p{qa)), and Ua ^{la-Qaa) Y{a<hKQab), 

respectively, which in turn will be abbreviated to f{q,0), {x'^s^\(j){q)) , and S{q — Q). 
Next, we can show[10] that the action of the total angular momentum J^p on 
can be expressed in terms of Cap on the collective amplitude f{qa,ds)' 



{xMi^} = J dT{x\6{q - QM{q))^apf{q, 0). (2.13) 
In this new representation, the operator Gap of (2.6) is rotated into Gab'- 



{AaS'ia\Gabma)) = \{l - ^„,) ^^^^^^^ ''^^^"^^ + 5aMaS'JGa\<t>{qa)) ■ 

^ qa-qb 



(2.14) 

Finally, using Eq. (2.11), (2.13) and (2.14), we have shown in Ref.[10] that the 
mean energy (2.4) is given in the body-fixed frame of reference by^ 

{m\i^)= J dTr{q,e){<f>iq)\Siq-Q)lH- KbMab] ' ^[ab]) 

a,b=l 

3 ^ . 1 . .... 

- Yl QaiGa - Pa) + 2 51 B[ab],[cd\{J[ab] " ^[ab]){J[cd\ " ■^[cd]) 

a=l abed 

+ X] B[ab],ciJ[ab] - ^[ab]){Gc - Pc) 
abc 

+ lY;^Ba,biGa-Pa)iGb-pb)]\m)fiQ,0), (2-15) 

^ ab ^ 



where 



^ _ Qab ^ 
Ha Hb 



iH, y[„b] , Qa = iH, Qa , (2.16) 



B[ab],[cdi = [Y[ab], iY[cd]] , B^ab],c = [Ylab], iQc] , Bab = [Qa^ ^<3b]- (2-17) 

Note that, in deriving the mean energy (2.15), we have terminated the series (2.4) at 
the quadratic terms in kap- This termination is justified by the validity of the adiabatic 
approximation in the present case, since we are dealing with nuclear dynamics for 



^ Recall that j[ab] and Ga operate on the intrinsic state \4>{qa)), but jC,[ab] and pa operate on 
the collective state f{qa,Os) (i.e., C[ab] acts on 6s and pa on qa). 
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which the time evolution of the collective variables is assumed to be slow on the scale 
of a single-particle (nucleonic) motion. 

As we are going to see next, the rotational and vibrational degrees of freedom 
appear explicitly in the energy expression (2.15); we will also show how to derive 
expressions for the collective rotational energy and for the moment of inertia. 

3 Description of Rotational States of Axially-Symmetric 
Nuclei 

Consider a permanently deformed, non spherical nucleus. Since we are interested in 
rotational motion only, we assume the nucleus to be in its vibrational ground state. 
In this case, we assume that the collective amplitude f{qa, ^) of (2.15) separates into 
a vibrational part, go{qa — Qa), and a rotational part, D{Cl): 

f{qa,n)=goiqa-qa)D{n). (3.18) 

The vibrational collective amplitude goiqa — Qa) represents here the zero-point oscil- 
lations about the equilibrium values, qa, of qa (a = 1, 2, 3). Hence, the wave function 
\'4)) of the system becomes {c.f. Eq. (2.3)): 

\^) = J drdn 5iq - Q)\<t>{qa))9o{qa - qa)D{Q) , (3.19) 

with dr = {qa - q2){q2 - g3)(^3 - qi)D^a=i'^qa a"d dO. = d^pdipsinOdO {c.f. Eq. 
(2.12)), and where 6{q-Q) is used to abbreviate Yll=i ^{qa - Qaa) Y{a<hKQab)- 
this case, after expanding \(t>{qa)) about \(j){qa)), the mean energy (2.15) becomes: 

{^\H\^) c^f dr I d^g*{qa - qa)D*{n){cj){qa)\S{q - Q)\h - ^ f „(4 - 4) 

J J ^ a=l 

+ \Y. B^b{Ja - Ca){Jb - Cb)]\'P{qa))g{qa - qa)D{n) 
^ ab ^ 

+ -^osc + -^coupb (3.20) 

where we have used the notation Oc to abbreviate {a, b, c being cyclic permutations 
of the body-fixed axes 1,2,3). In this expression, E^^^ is the energy of the zero- 
point oscillations; -Ecoupi is the coupling-energy between the rotational and vibrational 
motions which we can neglect. As for E^^^^, we will drop it from all following mean- 
energy expressions, since it represents only a constant shift of the entire (rotational) 
energy spectrum. Note that, in the derivation of the mean-energy (3.20), we have 
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ab\ 



ah) 



approximated the operator Bab by its mean value {4>\Bab\4') (i e., (^|-B, 
In what follows, the notation Ba will be used to abbreviate Baa 

For the sake of simplicity, we shall focus in ttiis work only on deformed, even, 
and axially-symmetry nuclei. Consider the axis 3, of the body-fixed frame, to be the 
axis of symmetry for the system. As a consequence of the axial symmetry, we have: 
9a = ^2 / qs and Bi = B2 = B B^. 

Now, since H, and Jz mutually commute^ commute, they possess joint eigen- 
functions. The structure of our trial function allows it to be an exact eigenfunction of 

and J2, but provides only a variational approximation to the energy. In the case 
of axial symmetry, this trial function |^) can be obtained from (3.19) by expanding 
\(j){qa))D{0.) in terms of the Wigner D— functions^: 



im) = J 



drdn S{q-Q)] 



121 + 1 
167r2 



\cl>K{qa))'D'MK{^) 



K 

9o{qa - qa) 



(3.21) 



+ {-lY\^^K{qa)-I^'M-K{^) 

where is an eigenfunction of J3 {i.e., J3\(j)K) = K\(j)K)) and Vj^j^ is an eigen- 
function to £3 {i.e., C^D^jj^ = KT>\jj^), C'^ and Cz- 't then follows that 

/V/m) =/(/+!) IV'/ Af) , Jz\^im) = M\iIjim) . (3.22) 

For the simpler case of the K = {) band, the wave function IV'/m) is given by: 



IM, 



-I 



drdn 5{q-Q)\ 



1 21 + 1 



87r2 



\Mqa))'i^Mo{^)9o{q - qa) 



(3.23) 



Note that (as a consequence of axial symmetry) the following important relation holds 
for both forms, (3.21) and (3.23), of \^)D{ny. 

{J3-C3MD{n)=0 . (3.24) 

In this case of axial symmetry, and after omitting E'^gj, + Ecoupi. we can see that 
the mean-energy (3.20) reduces to: 



(V'/mI^IV'/m) 



dTdng*iqa - qa)D*{Q){<t){qa)\5{q - Q) 

2 



B 



a=l 



\ct){qa))goD{n) 



+£1 + £2 



(3.25) 



* J is the total angular momentum and Jz is its Z component with respect to the Lab frame. 
®The definition of 'Dj^_i({^Y) used here is that of Bohr-Mottelson 



8 



with 

2 1 ^ 

H^ = H-Y^Yj^ + -BY,Jl , (3.26) 

a=l a=l 

r . 2 . _ 

Si = dTdng*D*{(t>\Siq-Q)Y,{Ya-BJa)Ca\(l>)go{qa-qa)Din) , 
■' a=l 

£2 = (S12 + B21) J dTdng*D*{cl>\S{q - Q){Ji - £i)(J2 - C2M)9oD{n) , 

(3.27) 

where \(j){qa))D{^) is given by (3.21) or (3.23), depending on whether one is inter- 
ested in the K ^0 band or the K = band. 

We should now specify the description of the intrinsic structure of the system. 
To this end, we assume that the intrinsic state |0(^a)) is given by a mean field 
approximation such that (<?!'(^a)|<3ii|<?!'(^a)) = qi is equal to (<^|(522|0) = q2 {is, 
such that qi = q2, the axial symmetry condition). This can be achieved by means of 
a constrained variational principle. 

Let us now look at the determination of the collective tensor operator Q. We 
determine the particle-hole (ph) components of the tensor operator Qab such that 

{4>\Ya — BJa\ct>) is variationally stable, i.e., 

5{ct>\fa-BJa\(t>)=^ (0 = 1,2) . (3.28) 

This variational condition insures that the simple expression (3.21) for 1-0) is adequate 

to describe the rotational energy (term ~ /Z^) correctly. 

To determine the mean energy (3.25), we need to calculate £1 and £2. In what 
follows, we are going to show that both £1 and £2 are identically zero. First, the term 
(^|(5(g — Q)(Ya — BJa)\4>) in the integrand of £1 can be rewritten as^ 

{mq-Q)&a-BJam ^ {mq-Qm{(t>\^a-B.ja\ct>) 

crfx 

(3.29) 

where a, r, . . . refer to unoccupied (particle) states, while i2,X,... refer to occupied 
(hole) states. Using the condition (3.28), we see that the term {(j)\S{q — Q){Ya — 

®In this approximate expression, we have neglected the two-body part of the operator Ya = 

[H,iYa]. 
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BJa)\(t>) becomes equal to {(j)\5{q — Q)\<j)){<j)\Y a — BJa\(l)) ■ Now, using the fact that 
\<j)k) is an eigenfunction to J3 and that the action of both Ya and Ja (a = 1,2) on 
\4>k) generate \(i)k±i) , we can ascertain that (</>|l"a ^ BJa\<j)) is itself identically zero, 
and hence £1 is equal to zero. To see this, note that {Ya — BJa) have non-zero 
matrix elements only between \4>k) and {4>k±i\- So, 

{(t>k\Ya-BJa\<t>k) = Q (a = 1,2) , (3.30) 

and also 

{(l)^k\Ya-BJa\(l)k) = ^ , (3.31) 

except for K = i but, in our case, K is always an integer (because we are dealing 
with an even nucleus). 

Second, £2 is identically zero, since both B12 and B21 can be shown to be equal 
to zero. To see this, using these expressions, 



Ji, iYi 



J2, iY2 



Q33 — Q22 

qs - q2 
Q33 - Qu 



Ji, iY2 
J2, iYi 



-Q12 

q2-qi 



(3.32) 
(3.33) 



93 - qi 

we can easily show the following important relation: 

(<^|[4, i%m=5ab (a = 1,2) , (3.34) 

since {Qab) = ^ab- Now, applying this relation to the variational principle (3.29), we 
can verify that (01 [Fa; ^feJI^) 's equal to BSab, i e., 

Bab = B{(l>\[Ja, i%m = B5ab (a = 1,2) . (3.35) 

Now, since both of £^1 and £2 are zero, and using the relation {J2,—tz)\(j))D{Q.) = 
of (3.24) and (3.28), we can show that the mean energy (3.25) reduces to 



H-^B{Jl + Jl + jl 



{^iMmiu) = J dTdng*o{nM6iq - Q) 

+1 1 drd^limq - Q)\<l>)9lD*{^) 



\4>)goD{n) 



Cl + Cl + Cl\goD(n) . 

(3.36) 



Using the approximation 



B 



By, 



{mq -Q)[H-^J']\<t>)^ {mq - Qm{4>\H - - J10), 



(3.37) 
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we can rewrite (3.36) in the following simpler form''^ 

= ~ri — Fi — r ~ o^^^^ + ^) + ^"^'^ ~ o^iJi +J2+ JsM) 



(3.38) 



where we have used he fact that C'^Vj^j^ = I{I + l)Vj^j^. Note that the energy 

expression (3.38) has a term, —\B{J'^) which represents a substraction of a mean- 
rotational energy. This term is familiar from the standard Peierls-Yoccoz angular 
momentum projection method. We expect this approximate treatment of the 5- 
function in (3.38) to overestimate the mean-energy by a term of the order of half the 
zero-point vibration energy. 

Moment of Inertia 

Let us now look at the moment of inertia, which is given by B"^. The inertia! 
parameter B can be determined from eqs. (3.29) and (3.34); i.e., it is given by the 
two equations 

{4>\ [h, iZ^ \(l>ph) = {<t>\Ja\<l^ph) , = (</,| [Za, iJa] |0) , (3.39) 

where Za = Ya/B. This expression for B^^ is of the well-known Thouless-Valantin 
form[13]. Note that if we neglect the residual two-body interactions from the Hamil- 
tonian, expressions (3.39) would give rise to Inglis cranking formula[14] 

1 . (3.40) 

This approximate formula is well-known to overestimate the moment of inertia quite 
badly. 

In what follows, we are going to examine the connection between the BO method 
and the large body of (rotational) data[15], e.g., the moment of interia increases 
with angular momentum /. First, note that Eq. (3.38), which was derived for time- 
reversal invariant |(/)), describe a rotational spectrum with constant moment of inertia, 
B~^ , in disagreement with data. To see this, consider the case = for which 
(^o|«^i|<!^o) = (i = 1,2,3). Hence, the energy expression (9.2.20) becomes 

Ei = ^BI{I + l) + {H)-^Efi , (3.41) 



'^Recall that we have omitted the vibrational energy part, E^^^ 
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with A£^/^ = y(^<^i + ^J^)' the fluctuation energy which is generated by angular 
momentum fluctuations. In this case, therefore, the energy spectrum is that of a 
rigid rotor, since the moment of inertia B~^, as given by (3.39), is constant. This 
contradicts, of course, the experimental facts. 

Second, note that the failure of (3.38) to generate a moment of inertia, B~^, 
which increases with angular momentum is due to a restrictive assumption on 
the time-reversal invariance of |^) for K = 0. In what follows, we are going to 
show that the BO approach has a natural mechanism for introducing a moment of 
inertia which increases with /, provided the restrictive assumption on \(p) is dropped. 
In addition, we will show that the energy we obtain for this case is lower than the 
energy, (3.41), obtained with a time-reversal invariant |^). To this end, let us consider 
a symmetry-violating for which ( Ji) is not zero but for which (J2) and (J3) are 
both zero. In this analysis, we will restrict ourselves to the simplest case: K = 0, 
and hence |(/)o) is an eigenfunction to J3 with eigenvalue zero, J3\4>o) = 0. In this 
case, the energy expression (3.38) reduces to 

Ei = ^BI{I + 1) + {H) - |(Ji)2 - AEfe . (3.42) 

Note that this energy is lower than the energy, (3.41), obtained with a time-reversal 
invariant \(J)) (provided the fluctuation energy is unchanged). In what follows, we 
shall neglect the angular momentum fluctuations, A.Ef£, in the determination of the 
mean field The energy expression (3.42) then provides a basis for a variational 
determination of the symmetry violating 10): 

S{H-ioJi)^ = , u; = B{Ji)^ , (3.43) 

where the notation (0)^^ is used to abbreviate (0(a;)|O|0(a;)). The parameter u has, 
obviously, the significance of an angular velocity. 

Now, we are in a position to show that B{lo) decreases when the angular momen- 
tum increases. To see this, using the relation B{io) = jfy- of (3.43), we have 

dB{io) _ 1 

Since {Ji)u) is well-known to increase with the angular velocity u, and as shown in 

Fig. 1, the slope ^^^^ is always larger than Thus, the slope, ^ of B{lo) is 

negative and, hence, B{uj) would behave as shown in Fig. 2: B{uj) decreases as the 
angular velocity 00 increases. 



1 - 



d{J: 



1 w 



{J,)/iO 



(3.44) 
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{Jl)u> - 




Figure 1: (Ji)a; increases with the angular velocity uj. 



Therefore, we conclude that the moment of inertia B~^{io) increases, indeed, with 
angular velocity lo, and hence with angular momentum also. 

In what follows, we are going to show that there exists a non-zero value, ujc, of 
oj at which the energy Ej{ojc) of (3.42) is equal to its lowest value. To this end, let 
us write the energy expression (3.42) in the following form (from which we omit the 
fluctuation term, ^Efi): 



Ei{lo) ~ {H-LoJi)^ + -B{uj)[{Ji)l + I{I + l) 



{H - ujJi)^ + - 



(3.45) 



First, note that the derivative, 

dEiicu) 1 



dB{uj) 



dio 



of Eiiuo) vanishes at a value ujc which is determined by (Ji)^^ = I{I + 1), i e 
dEi{oj 



(3.46) 



dw 



= 



(^1 



/(/+1) . 



Second, we can easily show that the second derivative of Ei{uj), 



(fEi{u) 
du^ 



UJ 



B^ico) 



LO f dB{uj) 



B{lo) \ dw 



dB{u) 
dco 



(3.47) 



(3.48) 
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Figure 2: The inertial parameter B{uj) decreases as the angular velocity u in- 
creases, since dB/dw < 0. 



is positive, since, as shown above [c.f. Eq. (3.44)), ^ is negative. Finally, we 
conclude that, using a trial function \4>{uj)) whose time-reversal symmetry is broken, 
one obtains, indeed, lower values for the inverse moment of inertia, B{uj), and for 
the energy than those calculated with a T ■ R invariant mean field. 

Calculation of the energy difference: AEj 

Let us now calculate the energy difference, AEj, between Ej{u = 0) and Ei{uc)- 
Using Eq. (3.46), we can show that 



AEi = Ei{loc) - Ei{0) = / 

Jo 



dEiito) , 

; du 

dio 



This expression can, after a partial integration, be reduced to 



(3.49) 



AEi = -^BoI{I + 1) + ^ ((^i)) '^('^^i) 



V^w) /, ^ ,x r B 



Jo 



B{u) 



d{Ji) . (3.50) 



This expression shows that Ei{u)c) is, indeed, lower than £^/(0), since Bq is larger 
than B{ijj). So, if we know the dependence of the angular velocity oo on {Ji)oj, we 
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Ei{u) + AEfi 




Figure 3: Behavior of Ei{u) as a function of the angular velocity where | A£^/| = 
\Ei{Q) — Ei{ujc)\ is an increasing function of the angular momentum I 

can easily calculate the energy difference between -5/(0) and Ej{iOc)- Note that, the 
energy difference |A£^/| is an increasing function of the angular momentunn /. The 
qualitative behavior of the energy Ej{ijj), for various values of/, is plotted in Fig. 3. 

4 Discussion and Conclusions 

Let us summarize what we have achieved in this work. First, we have shown that the 
moment of inertia generated by the NBO method is identical to the Thouless-Valantin 
form. Second, the two relations (3.43) and (3.47) determine the intrinsic (symmetry- 
breaking) function |(/>t^) and the value, lOc, of uj where Ei{ujc) is equal to its lowest 
value, respectively. These two relations provide a bridge (connection) between the 
NBO method, which is a truly quantum mechanical description of collective motion, 
and the semi-classical approaches based on the idea of self-consistent cranking (SCC). 
Thus, we have established a connection between the NBO approach and the large 
body of experimental data, since the two relations (3.43) and (3.47) are known to 
provide reasonable descriptions of vast amounts of empirical data ranging from low- 
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lying rotational states to high angular momentum states. [16]-[21] So, the present 
(NBO) method appears to be well-equipped to describe low as well as high lying 
rotational states. Additionally, we should mention that work has been started to 
apply the NBO method to the description of the backbending phenomenon which 
was first observed by Johnson and his collaborators. [22] 

In summary, we have studied here the rotational spectrum of even, axially-symmetric 
nuclei within the framework of the NBO method. We have made use of trial functions 
in which the intrinsic structure is described within a mean-field approximation. We 
have shown that the NBO formalism gives back the Thouless-Valantin moment of 
inertia. Then, we have established a connection between the NBO method and the 
sec model, which has been successful in reproducing vast amounts of experimental 
data. Finally, we have shown that the introduction of a non time-reversal invariant 
intrinsic function both lowers the energy for a given /, and provides a moment of 
inertia that increases with the angular momentum I. 
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